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We present a method to derive explicit forms of tight correlation function Bell inequalities for 
three systems and dichotomic observables, which involve three settings for each observer. We also 
give sufficient and necessary conditions for quantum predictions to satisfy the new inequalities. 
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(3 ' INTRODUCTION 
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The problem of the possibility of a local realistic interpretation of quantum mechanics was first addressed in the 
^ ' discussion between Einstein, Podolsky, Rosen [l| and Bohr 01 • The final answer was provided in 1964 by Bell 0|. 
The paper of Bell contains an inequality for local realistic (LR) correlation functions, which can be violated only by 
non-classical, entangled states. The correlation function is defined as a mean value of a product of local observables. 
] For example, a two-qubit correlation function shall be given by Eij — {AiBj), with Ai being an zth observable for the 
first observer (AHce) and Bj is jth for the second (Bob). The observables are dichotomic, with eigenvalues ±1. In the 
' work of Clauser, Horne, Shimony, and Holt (CHSH) Q a new inequality was derived, which, in contradistinction to the 
' original Bell's expression, can be applied to real physical processes. 21| Again, it was formulated for LR correlation 
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, functions. Such formulations of the Bell inequalities are now the most popular forms of it. The discovery of the 
Greenberger-Horne-ZeiHnger correlations triggered generalizations of Bell inequalities to more subsystems [5y. The 
early efforts in this direction were restricted to the scenarios involving two experimental settings per observer only. 
The full set of such inequalities was finally given in Refs. (1, 01- More recently, however, maps between general Bell 
inequalities for correlation functions (even yet unknown ones) and quantum cryptography as well as quantum 
communication complexity problems [9||, have been found. Thus, there is an additional motivation to find new tight 
[ Bell inequalities of this kind, with the obvious direction of generalization toward more allowed settings for each 
observer. Allowing more settings leads e.g. to more complicated communication complexity problems, and from the 
fundamental point of view, must further restrict the realm of classically admissible correlations. 

■ By tight Bell inequalities one understands these, which define walls of a polytope in a certain statistical hyperspace, 
^ ' which in turn defines the full set of possible local-realistic predictions for the values of correlation functions for the 

given set of Bell-type experiments. The vertices of the polytope are deterministic predictions, that is, are given by 
J> ' extremal correlation functions for which the product of the values for the each pair of measurement directions are 
k>( always the same. Thus vertices have components equal to ±1. The polytope was initially proposed for the probability 
rS , hypers pace lldl . Various aspects of the problem were subsequently discussed by a number of authors, for example 

■ in [ill , lid . Il3l |. More recently, a method to derive the general structure of coefficients of tight Bell inequalities for 
correlation functions, which involve three, or more, settings available to each observer, was presented in [l3|. Here, 



basing on the results of we develop this method with the aim of getting the explicit form of the inequalities for 



three observers. This is the simplest case for which tight three-setting Bell inequalities are thus far unknown. 



THREE QUBIT INEQUALITIES 

Let us consider a situation, in which each of the three observers can choose between three measurements. For three 
systems (observed by Alice, Bob and Carol independently) and dichotomic observables, the tight Bell inequalities 
have a general form of [14 1 

- 1 < ^ J2ai,b,,ck=±i '5'(ao, ai, 02; &o, 6i, &2; Cq, Ci, C2) 

X (EL,fc=o««^jCfe£;,jfe) < 1- (1) 

where Eijk is the value of the local realistic correlation function (i.e. the average of the product of the local results, 
which in turn are assumed to be ±1) for settings k of A, B, and C. S{...) is a sign function (values ±1), of the 
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nine parameters, oq, ...,C2. An admissible sign function has the form of [lj| 

S{ao,ai,a2; foo, ^i, 62; cq, ci, C2) 
= EL,fc=o 9tjkarb-jCk = ±1. (2) 

This relation can be inverted to give 

9ijk = ^ S{ao,ai,a2;bo,bi,b2;co,ci,C2)aibjCk. (3) 

This allows rewriting ([1]) as: 

2 

X! < 1, (4) 

i,j,k=0 

where iJy/c = (AiBjCk). 

Note also that the form of the sign functions, as in equation restricts the admissible class which generates the 
inequalities. The elementary trait of the class is that the admissible sign functions do not depend upon product of 
indices pertaining to one observer, e.g. aia2, etc. Note also, that one can reduce the number of considered inequalities 
to those involving sign functions with qq — 1, and in such a case the normalization factor changes from ^ to This 
convention was used in ref. |:14]. However, it turns out that the more symmetric convention is handier for pinpointing 
the exact form of the coefficients, and thus we shall adopt it here. 

PROPERTIES OF THE COEFFICIENTS IN THE INEQUALITIES 

To derive explicit Bell expressions it is convenient to define the following quantities, which we will further on call 
deltas. The first order deltas: 

A, = l(5(x = l)-5(x = -l)), (5) 

the second order deltas: 

A,^ = i(A,(y = l)-A,(2; = -l)), (6) 

A,y, = ^(A,y(z = l)-A,y(z = -l)). (7) 

Each of the indices y, z is related to a different party, e. g. a; = ao, ai, 02; y = h^, 61, 62; z = cq, Ci, C2. 

The third order deltas are equal to the coefficients gju- We can create the second order deltas from them in the 
following way: 

(co,Ci,C2) = Co 

and likewise for the remaining l^aibj 's. Likewise, 

Aao(^0, ^1, ^2; Co, Ci, C2) = 6oAaofco(cO, Cl, C2) + 6lAao&i(co, Ci, C2) + h2^aob2{cQ, Ci, C2), 

(9) 

and 

<5'(ao, oi, 02; 60, ^1, ^2; Co, Cl, C2) 
= aoAao(6o,&i,&2;co,ci,C2) + aiAai(6o,6i,62;co,ci,C2) +a2AQ2(6o,foi,62;co,ci,C2) , (10) 



and the third order deltas: 



etc. 
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Properties of the deltas 



As 5 = ±1, each of the first order deltas Aa^ can be equal to ±1 {S changes with a^) or (5* does not change with 
Qi). From the possible values of the first order deltas we can deduce, that a second order delta can only be (no 
change), ±^ (a change of Aq. with from to ±1, or vice versa), or ±1 (a change between 1 and —1). By the same 
argument the coefficients gijk may assume values 0, ±| or ±1. 

Finally, one can easily see that the coefiicients gijk satisfy the following relation 



E 

j,k,l=Q 



9]kl 



The first condition is due to a fact that the inequalities we are searching for are saturated for all deterministic LR 
models, the second is related to Parseval's theorem: for every sign function of nine dichotomic variables we have 

The last inequality immediately implies that when one of coefficients gijk has modulus 1, then all the others must 
be 0. In such a case, we get a family of trivial inequalities \Eijk \ < 1. 



Construction of sign functions 

Each second order delta can contain one, two, or three coefficients gijk, and, similarly, the first order deltas contain 
up to three second order deltas. Not every pair or triplet of the second order deltas can be arranged into higher-order 
deltas, however. The following table lists the families of Aa^bj 's, and what other Aa^tj 's they can go with to create a 
valid Aai (we avoid repetitions in the table): 



Aoifij 


goes with: 


i((-l)"3c, + (-!)%) 


±i((-l)"c, -(-!)%) 


i(2(-l)™c, + (-l)% + (l)°c.) 


±i((-l)™c, - i-ircy) or ±i(2(-l)™c, - (-l)"c, - (-l)°c,) or 
±i((-l)™c, - i-ircy) and ±i((-l)™c, - (-l)°c,) 


i((-l)"c, + (-l)"c,) 


alone or with ^{{—l)"^Cx - {-l)^Cy) or 
two i((-l)"c. - (-l)"c„)'s 




any one of ±^Cx, ^\cy, ±\cz or ±j{{—l)"^Cx + (-l)"cy) and 
±i((-l)"c, - {-iTcy) or ±\{{-lTcy + (-l)°c,) and ±K(-l)"Cy " i'^Tc.) 


i((-l)™c, + (-l)"c,) 


±i((-l)"c, + {-iTcy) or ±\{{-\Ycy - (-l)°c,) and ±i((-l)'"c, + (-l)°c,) 



In the table {x, y, z) stand for a permutations of (0, 1, 2), while 711,71,0 ^ {0, 1}. 

Thus any first order delta belongs (after some transformations, i. e. permutations of observables, or sign flips) to 
one of the families listed below: 





Aa. 


|Aa.P 




Aa. 


|Aa.P 


Ao 


5(&o(co + Cl) + 61 (co - Ci)) 


16 
Ifi 


Ai 


i(6o(-3co + Cl) + 61 (co + Cl)) 


12 
Ifi 


A/7 


i&o(co + Cl) 


8 

Ifi 


Am 


i(6o + 6i)co 


8 

Ifi 


Aiv 


i(6oco + 61C1) 


8 
Ifi 


Av 


i(6o(2co + Cl + C2) + 61 (ci - C2)) 


8 

16 


Avi 


i(feo + ^'i)(co + ci) 


4 
Ifi 


Avii 


i(6o(ci - C2) + 5i(co - Cl) + 62(00 - C2)) 


6 
16 


Aviii 


i(&o(2co + Cl + C2) + 61 (co - Cl) + 62(00 - C2)) 


10 
Ifi 


Aix 


i(6o(2co + Cl + C2) + 6i(2co - ci - C2)) 


12 
16 


Ax 


i(5o(2co + 2ci) + (61 + 62) (co - Cl)) 


12 
Ifi 


Axi 


i(26o + &i(co + Cl) + 62(00 - Cl) 


8 
16 


Axil 


i(26oCo + 61 (ci + C2) + 62 (ci - C2)) 


8 
Ifi 









The norm |Aa. p is the sum of squares of those gijk coefficients, which enter a given delta. As it has been already 
mentioned before, in an admissible sign function these norms add up to 1 jloj. However, this is only a necessary 
condition to build an inequality. The necessary and sufficient condition is that X)i=o lAajP = 1 and Aa^Aa., = for 
i ^ i'. Only a few pairs or triplets of these deltas lead us to a proper Bell expression. 
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TWO QUBITS 

As an illustration of the method we will consider the two qubit case first. 

For two systems (observed by Alice and Bob), the tight Bell inequalities have a general form of [iJ] 
-1< ^Ea.,&,=±i'S'(ao,ai,a2;5o,5i,52) (EL=o«i^J^y) ^ ^■ 

(11) 

respei 

S{ao,ai,a2;bo,bi,b2) = Z]i,j=o = 

(12) 



where Eij — (AiBj) is the value of the correlation function for the ith and jth setting of Alice and Bob, respectively. 
S{...) is the sign function of ao, ai, 02, 6o, 61, 62- The coefficients 5^ 's of this function are related to S by 



Thus 



This allows to rewrite as 



9v ^ X! S{ao,ai,a2;bo,bi,b2)aibj. (13) 

ai — ±1 



i,j=0 



< 1. (14) 



It is enough to define first and second order deltas. The second order deltas are equivalent to gij's. As we have 
mentioned, the first order deltas can only take values of or ±1 and g^j — 0, ±i, or ±1. If some Aq^ is always equal 
to ±1, the sign function is factorisable, and generates a trivial bound: 

|Sy|<l. (15) 

If some takes all three values (±1 and 0), it can be locally transformed to ^(^0 + ^i)- The other A^. vanishes 
when the first one does not, and vice versa, so it can be taken as ^{bo — &i). Thus we reproduce the CHSH inequality 

i 1^00 + ^^01 +^10 -£^11 1 <1- (16) 

No other inequalities of this kind for two qubits are allowed, because we have utilized all possible values of g^-'s and 
all possible forms of Aq^ 's. 

As we see the set of Bell inequalities for correlation functions, for two partner and three settings on each side, boils 
down to the good old CHSH inequalities, for all possible combinations of pairs of settings. This result corroborates 
with the earlier finding of Garg [iSj. 

THREE-QUBIT INEQUALITIES 

Since the method is universal, it produces both the well-known standard Bell inequalities 0,0], which are shown 
in the first subsection below, and the new ones, which are the subject of the second subsection. 

2x2x2 inequalities 

• Using A/ and A/y we can obtain the following sign function: 

i(ao(5o(-3co + ci) + 6i(co + ci)) + ai(6o + fei)(co + ci)), (17) 
and by putting it into ([3]), we obtain: 

-| — 3i?ooo + £-001 + £-010 + £-011 + £-100 + £-101 + £110 + £111 1 < 1- (18) 



5 



Taking two A/7's or two A///'s we get 

i(ao6o(co + ci) + ai^o(co - ci))- (19) 
This leads to the inequality, which is a trivial extension of the CHSH one: 

21-^000 + £'001 + -Eioo — £-101! < 1, (20) 
With two A//'s or two A///'s one can also get 

i(ao(6o + 6i)co + ai(6o - &i)co), (21) 



this gives 



• and 



this in turn implies 



21-E'OOO + Eoio + Eioo — Eiio\ < 1. (22) 



"S* ^(ao&o(co + ci) + ai&i(co - ci)), (23) 



21-^000 + £'001 + -Eiio — I < 1. (24) 



Finally, a combination of A// and A/// leads to 



S = i(ao(6oCo + ^ici) + ai(6iCo - 60C1)). (25) 
It results in a Mermin-Ardehali-Belinskii-Klyshko inequality [l^ : 

-\Eooo + Eon + Eioi — Eno] < 1. (26) 

This closes the group of possible inequalities which involve only two settings (out of the three ones allowed in the 
considered scenario). 

Inequalities with up to three measurement settings per observer 

However, our method can also lead to thus far unknown inequalities. The most important is one in which all three 
observers perform three alternative measurements. The sign function is built of Ay and two Ay//'s: 

S = j{ao{bi{ci - C2) + 62(ci - C2)) + ai(6o(ci + C2) + &i(co + ci) + 62(-co + C2)) 

+a2(&o(-ci - C2) + &i(co + C2) + 62(-co + ci))). (27) 

It represents the following inequality: 

i|£'oii — -E012 + £'021 — Eq22 + Eioi + E102 + EiiQ + Em 

~Ei20 + -E/122 ^ -£'201 ~ -£'202 + -£'210 + -£'211 ~ -£'220 + -£'221 1 !i 1- (28) 

From this inequality follow two other, less general conditions. Putting oq = 02 we get 

S = i(-ao(6o(ci + C2) - 6i(co + d) + &2(co - 2ci + c2) 

+ai(&o(ci + C2) + 5i(co + ci) + 62(-co + C2))), (29) 

this represents 

■J I ~ -£-001 ^ -^002 + -EqIO + Eqii — -B020 + 2i?021 

+-E022 + EiQi + E102 + EiiQ + Eiii — E120 + -B222I < 1- (30) 
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If we rather put 02 = ai we have 



this gives 



S = |(ao(6i(ci - C2) + &2(ci - C2) 
+ai(6i(2co + ci + C2) + 62(2co - ci - C2))), 



zl^-Oll ^ £'012 

£■112 - 



4 

+£^111 



- £'021 — -£022 + 2i?iio 

2i?120 — -£121 — £'122 I < 1- 



Our method leads to yet another Bell inequality. Two Ax//'s can be arranged to 

S = i(ao(26oCo + 61 (ci + C2) + 62(21 - C2)) 
+ai(26oCo - 61 (ci + C2) - 62(01 - C2))), 



what represents 



3:|2-Booo + £-011 + £-012 
+2£'ioo — £'111 ~ £-112 



- Eq 

E121 



- Eq22 

£^122! < 1- 



Interestingly, inequalities l(30|) , l(32|) , and l(34|) are special cases of the one given in [l7|, \l 



g:|£'Ooo 

+Eq22 - 

+E100 - 



- Eqoi 

E023 
Eloi 



-£'122 ~ £^1 



23 



" £^010 — Equ 
£^032 — £'033 

£-110 — Em 

£'132 + £'133 < 1- 



(31) 



(32) 



(33) 



(34) 



(35) 



Inequality (|34l) is obtained by putting Cq — Ci. Then we can get l|32p by ]3q = B2. Finally, choosing Bo — B2 and 
Cq — C2 from ((35l) we reach ((30|) Further such simplifications lead to the standard 2x2x2 inequalities. 



CONDITIONS ON STATE TO SATISFY THE NEW INEQUALITIES 



In this Section we will derive sufficient conditions which a state must satisfy in order to never violate the Bell 
inequality. Such conditions are not dependent on particular measurement settings, but refer only to the properties of 
the state. 

As inequalities (|30l) and l(32|) are special cases of l(28l) . we focus only on the conditions for i(28|l and (|34l) . 
Let us recall that any observable on a qubit with spectrum { — 1, 1} can be written as Xi = Xi ■ {X = A, B,C; 
i = 0,1, 2, the upper index of the Pauli matrix vector enumerates the qubit and all the vectors are normalized). In 



terms of the correlation tensor (elements of which are given by Tijk 
written as 



Trp{af ® erf (g) )) inequality l|28p can be 



C2 



+ {Ai - A2) ® -Bo ® {Ci + C2) 
+ {Ai + A2) ® (Bi + B2) ® Co 
-\{Ai + A2) ® {Bi + B2) ® (Ci + C2) 

i(Ai - A2) ® (Bi - §2) ® (Ci - C2)) 



< 1, 



(36) 



where the tensor scalar product is taken as ^ i? = j k=i ^ijkBijk Note that the tensor in (|36|l constructed from 
the vectors defining the measurements contains five mutually orthogonal (in the sense of 0) terms. 
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We now choose local coordinate systems, in which 

Ai + A2 = 2cosaei, (37) 

Ai - A2 = 2smae2, (38) 

Bi +^2 = 2cos/3ei, (39) 

Bi - B2 = 2sin/3e2, (40) 

Ci +(72 = 2 cos 761, (41) 

Ci- (72 = 2 sin 762, (42) 

where ei = (1,0,0), 62 — (0,1,0), and 63 — (0,0,1). Let us moreover introduce a short-hand notation Taij = 
T QAo^Ci^Cj, TiBj = T Q Bi (g) Bq ^ Cj , and Tijc = T (D Ci® Cj ® Cq. This allows to rewrite (|36|) in a form of scalar 
product of two real vectors: 

|(7ai2, ?2Bi, Ti2C, Jill, T222) 

■(cos/?sin7,sinacos7,cosQ;sin/3,cosacos/3cos7,sinasin/3sin7)| < 1. (43) 

This is the necessary and sufHcient condition for the given state to satisfy the inequality. That is, if the maximum of 
the left hand side for all possible local coordinate systems, vectors Aq, Bq and Co, and angles a, (3 and 7 is less or 
equal to 1, then the predictions for given state always satisfy the inequality. 

If one employs the Cauchy inequality, \vi ■ V2\'^ < |wi 1^1^21^, one can formulate a concise sufficient condition on a 
state to satisfy (|28l) . The state cannot violate the inequality if 

Thi + + Tl,2 + Tiei + T?2C < 1 (44) 

holds in all local bases and for all Aq, Bo,Cq. Finally, we utilize the fact that = X]i=i ^^12 (this is 

because the optimal vector is Aq^* = (rii2, T212, T312)/ \/Tii2 + T212 + T^i2)j etc., and obtain the following form: 

3 

Tfn + TI22 + E(^'i2 + + Tl2^) < 1. (45) 

i=l 

Let us derive a similar condition for (|34| . In the fashion of l(43l) the inequality reads 

|f (cosaei (g) So «> (7o 
+ sin Q!e2 <Xi (cos /3 cos 76*1 (g) e 1 + cos (3 sin 76*1 (g) 62 

+ sin/3cos7e2 (8) ei — sin/3sin7e2 (g) 62)1 <1 (46) 

(we have renumbered observables of Alice). Now, repeating the above argument, we get a necessary and sufHcient 
condition on a state to satisfy l(34|) . Namely in all local coordinate systems and for all Bq and Co we shall have 

2 

T?BC + E ^2,, < 1, (47) 



or using the trick presented earlier 
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CONCLUSIONS 



A method to derive explicit Bell inequalities with more then two observables per site was proposed and demonstrated 
for the case of three particles and up to three observables per site. New tight Bell inequalities for correlation functions 
and the respective conditions to satisfy then are found. In particular, we present an inequality, ((28)) . in which all three 
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observers perform three alternative measurements. Interestingly, the inequalities l|30p . (|32|) . l(34|) can be shown to be 
implied by both the ones found in 17, l3| and (|28| . 

The method presented here can be generalized to M > 3 settings Bell inequalities, and more than three parties. 
This definitely would lead to narrowing the class of states which allow a local reahstic model. An old result of one 
of us [31, which shows a drastic constraint on such a class if the number of settings goes to infinity (for 4 parties or 
more), clearly supports this conjecture. 
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